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The disruption of the NRF2 (nuclear factor erythroid-derived 2-like 2)/glutathione-mediated antioxidant defense 
pathway is a critical step in the pathogenesis of several chronic pulmonary diseases and cancer. While the 
mechanism of NRF2 activation upon oxidative stress has been widely investigated, little is known about the 
endogenous signals that regulate the NRF2 pathway in lung physiology and pathology. Here we show that an E-box- 
mediated circadian rhythm of NRF2 protein is essential in regulating the rhythmic expression of antioxidant genes 
involved in glutathione redox homeostasis in the mouse lung. Using an in vivo bleomycin-induced lung fibrosis 
model, we reveal a clock "gated" pulmonary response to oxidative injury, with a more severe fibrotic effect when 
bleomycin was applied at a circadian nadir in NRF2 levels. Timed administration of sulforaphane, an NRF2 
activator, significantly blocked this phenotype. Moreover, in the lungs of the arrhythmic Clock A19 mice, the levels 
of NRF2 and the reduced glutathione are constitutively low, associated with increased protein oxidative damage and 
a spontaneous fibrotic-like pulmonary phenotype. Our findings reveal a pivotal role for the circadian control of the 
NRF2/glutathione pathway in combating oxidative/fibrotic lung damage, which might prompt new chronother- 
apeutic strategies for the treatment of human lung diseases, including idiopathic pulmonary fibrosis. 
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In mammals and humans, the circadian clocks in the 
brain and periphery are responsible for generating ~24-h 
daily rhythms that control physiology, metabolism, and 
behavior (Reppert and Weaver 2002; Hastings et al. 2003). 
At the molecular level, circadian rhythms emerge from 
the interlocked transcriptional/translational feedback 
loops consisting of core clock genes and oscillatory 
metabolic products (Ko and Takahashi 2006). The core 
clock transcription factors CLOCK and BMAL1 activate 
the expression of repressors Period (Per) and Cryptochrome 
(Cry) genes as well as a spectrum of clock-controlled genes 
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(CCGs). Subsequently, accumulation of PER/CRY protein 
complexes inhibits the BMAL1 /CLOCK activity, leading 
to rhythmic repression of their own transcription and 
of other CCGs. Many of the known CCGs encode rate- 
limiting enzymes and transcription factors that play 
pivotal roles in cellular metabolism and function. None- 
theless, the battery of genes under circadian regulation is 
largely tissue-specific so that circadian clocks modulate 
physiological processes unique to each organ (Panda et al. 
2002; Gossan et al. 2013). 

Daily changes in the environment (cycles in light/dark- 
ness, feeding, rest-activity, and temperature fluctuations) 
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inevitably expose the mammalian cells and tissues to 
periodic challenges, including oxidative insults from en- 
vironmental toxins/pollutants and endogenously pro- 
duced reactive metabolites as products of respiration 
(Patel et al. 2014). Failure to scavenge reactive oxygen 
species (ROS), reactive nitrogen species (RNS), and re- 
active electrophiles results in excessive oxidative stress, 
leading to damage to critical macromolecules and cellu- 
lar structures (Veal et al. 2007). The ability to anticipate 
and withstand such cyclic insults (e.g., through effective 
ROS scavenging) is therefore essential for normal pro- 
tective tissue functions. There is emerging evidence that 
circadian clocks regulate processes that keep ROS at 
physiological levels and protect organisms from oxida- 
tive stress (Stangherlin and Reddy 2013; Patel et al. 
2014). Indeed, self -sustained rhythms in the cellular 
redox state (Kondratov et al. 2006; O'Neil et al. 2011; 
Edgar et al. 2012; Wang et al. 2012) as well as ROS- 
responsive antioxidant genes (Geyfman et al. 2012; Lai 
et al. 2012) have been observed in the brain and/or 
peripheral tissues. Moreover, diurnal rhythms in or- 
ganismal survival to ROS-generating agents have been 
reported (Krishnan et al. 2008). However, the underlying 
molecular mechanisms of clock-controlled regulation 
of protective antioxidant responses and their role in 
tissue-specific physiology and pathology are still not well 
understood. 

One of the major cellular antioxidant defense pathways 
is mediated through the cap'n'collar basic leucine zipper 
transcription factor NRF2 (nuclear factor erythroid-de- 
rived 2 -like 2). The current paradigm of NRF2 regulation 
has mainly focused on oxidative stress-induced post- 
translational mechanisms, including regulation by its 
negative regulators (KEAP1, BACH1, and (3-TrCP), posi- 
tive regulators (DJ1, p62, and p21), or protein modifica- 
tions (Itoh et al. 1999; McMahon et al. 2003; Chen et al. 
2009; Komatsu et al. 2010; Chowdhry et al. 2013). In 
response to oxidative stress, NRF2 protein translocates to 
the nucleus and binds to the antioxidant response ele- 
ments (AREs) in the promoters of many antioxidant 
genes. Subsequently, NRF2 induces a transcriptional pro- 
gram that maintains cellular redox balance and protects 
cells from oxidative insults (Reddy et al. 2007; Malhotra 
et al. 2010). NRF2 transcriptional targets include gluta- 
thione cysteine ligase (GCL), involved in the rate-limiting 
step of biosynthesis of antioxidant glutathione,- glutathi- 
one-utilizing enzymes such as glutathione S -transferase 
(GST); and haeme oxigenase 1 (HMOX1), involved in the 
catabolism of pro-oxidant haeme. 

Nr/2-deficient mice show a decrease in both constitu- 
tive and inducible expression of multiple glutathione- 
dependent enzymes (McMahon et al. 2001), associated 
with an impaired antioxidant defense and susceptibility 
to stress-induced tissue pathologies (Cho et al. 2006; Chan 
and Kan 1999). Consequently, they are highly sensitive to 
many environmental toxicants and chemical inducers of 
oxidant injury and fibrosis in several tissues, including 
bleomycin-induced and butylated hydroxytoluene (BHT)- 
induced pulmonary fibrosis (Chan and Kan 1999; Cho et al. 
2004; Kikuchi et al. 2010). In humans, altered NRF2 



expression and impaired redox balance have been 
associated with the pathogenesis of chronic lung 
diseases (asthma, COPD, and idiopathic pulmonary 
fibrosis) as well as lung cancer (Cho et al. 2006; Hayes 
and McMahon 2009). Indeed, pharmacological target- 
ing of NRF2 has emerged as a novel therapeutic ap- 
proach to combat oxidative stress seen in chronic lung 
diseases (Malhotra et al. 2011; Artaud-Macari et al. 
2013). 

Recently, mRNA levels of Nrf2 as well as several of its 
target genes have been reported to show diurnal variation 
in mouse livers and Drosophila heads (Beaver et al. 2012; 
Xu et al. 2012). However, whether and how endogenous 
circadian signals control NRF2-dependent antioxidant 
activity in tissue physiology and upon pathological 
challenges, especially in the lung, has not been defined. 
Here we investigated the circadian rhythms in the pro- 
tein levels of NRF2 and activity of the NRF2/GSH 
pathway in lung tissues from light/dark cycle-entrained 
mice. We also report the physiological and pathological 
importance for the circadian control of the NRF2/GSH 
pathway using a well-established lung fibrosis model 
following a bleomycin challenge. Our data reveal a piv- 
otal role for endogenous NRF2 rhythms in the lung in 
coupling protective antioxidant responses to the time- 
of-day susceptibility to oxidant injury and pulmonary 
fibrosis. 

Results 

Circadian protein expression of NRF2 transcription 
factor in the mouse lung 

To investigate the hypothesis that the NRF2-mediated 
antioxidant pathway is under circadian regulation in 
the lung, we first examined the total NRF2 protein 
levels in mouse lungs harvested at 4-h intervals from 
mice kept in constant darkness (DD). Western blotting 
using a specific NRF2 antibody (confirmed in Nrf2 
knockout mouse embryonic fibroblasts [MEFs]) (data 
not shown) revealed a robust rhythmic pattern of 
NRF2 protein (one-way ANOVA, Tukey test, P < 0.05). 
The NRF2 protein peaked at the circadian time 3-7 
(CT3-CT7), reaching a trough at CT15-CT19 (Fig. 1A). 
Nuclear accumulation of NRF2 protein is essential for 
its subsequent activation of downstream genes (Reddy 
et al. 2007; Nguyen et al. 2009). Indeed, we detected 
rhythmic levels of NRF2 protein in the nuclear lung 
fractions, with the maximal levels corresponding to the 
peak of the total pool of NRF2 protein (P < 0.05) (Fig. IB). 
Immunostaining of NRF2 revealed a significant time-of- 
day difference in the lungs that was mainly localized to 
the bronchial and alveolar epithelium, with an approx- 
imately threefold weaker signal at "zeitgeber" time 12 
(ZT12; lights off) compared with ZT0 (lights on) (Fig. 
1C). The rhythmic NRF2 protein expression was also 
evident in the total cell lysates (Fig. ID) and nuclei (Fig. 
1E,F) of clock-synchronized Rati fibroblasts, demon- 
strating the cell-autonomous nature of the NRF2 protein 
rhythms. 
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Figure 1. Circadian expression of NRF2 protein in mouse lungs and fibroblast cells. [A, B) Representative Western blots of rhythmic 
NRF2 levels in total lung lysates and nuclear extracts from wild-type (WT) mice kept in DD. NRF2 densitometry (mean ± SEM) was 
normalized to p-actin or histone HI, and the highest expression set as 100%. (Gray bar) Subjective day; (black bar) subjective night. 
One-way ANOVA for the effect of time, (*) P < 0.05. (C) Representative micrographs of histological lung sections showing temporal 
NRF2 immunoreactivity at ZTO versus ZT12. Brown NRF2 staining was localized mainly to bronchial epithelium (BE) and alveolar 
epithelium (AE). Mean ± SEM; (***) P < 0.001, t-test. Bar, 100 |jlM. (D) Representative Western blots of temporal NRF2 levels in total 
lysates of Rati fibroblasts following serum shock. NRF2 densitometry (mean ± SEM) was normalized to tubulin. Harvest of cells 
started at 12 h after serum shock to exclude the initial gene induction. (**) P < 0.01, one-way ANOVA for the effect of time. (£, F) 
Representative immunofluorescent micrographs and quantification of endogenous NRF2 in Rati cells at 24 h or 36 h after serum shock. 
Nuclei were counterstained with DAPI. Bar, 50 |xM. (***) P < 0.001, £-test. 



The circadian clock exerts transcriptional control 
of the Nrf2 gene via a putative E-box element 
in its promoter 

The circadian clock transcription factors regulate target 
genes through the cis-acting elements in their gene pro- 
moters (Panda et al. 2002; Rey et al. 2011). We therefore 
analyzed the proximal promoter regions of the mouse, 
rat, and human Nrf2 genes for clock-regulated transcrip- 
tional elements, which revealed a conserved putative 
E-box (CACGTG) close to their transcriptional start sites 
(TSSs; —704 base pairs [bp] relative to the TSS in the 
mouse Nrf2) (Fig. 2A). Next, we cloned the wild- type 
mouse Nrf2 ::luc promoter (1.7 kb) and generated an 
E-box mutated construct. Coexpression of circadian tran- 
scription activators BMAL1 and CLOCK in Rati cells led 
to a significant stimulation of Nrf2 : :luc activity (2. 1 -fold; 



P < 0.05), similar to the activation of the Perl ::luc 
promoter (containing three canonical E-box elements,- 
2.4-fold; P < 0.01) (Fig. 2B). Mutation of the core E-box 
sequence in the Nrf2 promoter completely abolished its 
induction by CLOCK/BMAL1 complexes (Fig. 2C). More- 
over, coexpression of circadian transcription repressors, 
either PER2 or CRY1, suppressed the CLOCK/BMAL1 
transcriptional activation (Fig. 2D). Finally, stable trans- 
fection of the Nrf2 ::luc reporter into Rati cells revealed 
a clearly discernible circadian rhythm with a period of 
23.7 h (data not shown). 

We next examined the temporal profiles of endogenous 
Nrf2 mRNA in the lungs of wild-type and Clock A19 mice 
(a model in which the positive arm of the clock is ge- 
netically disrupted) (Vitaterna et al. 1994). This revealed 
a clear rhythmic expression of Nrf2 transcripts in wild- 
type lungs (P < 0.05, one way ANOVA) (Fig. 2E). In 
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Figure 2. Orcadian clock components control Nrf2 transcription via a putative E-box element in vitro and in vivo. (A) A conserved 
putative E-box element (shown in red) in the proximal promoters of the mouse, rat, and human Nrf2 gene. Sequences were aligned using 
the Clustal W algorithm and shaded with the Boxshade algorithm to highlight the conserved regions (shown in black). [B, C) Effects of 
CLOCK/BMAL1 overexpression on wild- type (WT) or AE-box iVr/2 ::luc promoter activity in Rati cells. Data (mean ± SEM) were 
normalized to the (3-galactosidase reporter gene and expressed relative to pCDNA3 control. The Perl ::luc reporter served as a positive 
control. (*) P < 0.05; (**) P < 0.01; (ns) not significant, t-test. [D] Effects of PER2 or CRY1 overexpression in suppressing the CLOCK/ 
BMAL1 -activated Nrf2 ::luc promoter activity in Rati cells. (*) P < 0.05; (**) P < 0.01, t-test. [E] Temporal mRNA expression of Nif2 in 
wild-type and Clock A19 lungs. Data (mean ± SEM) were normalized to Gapdh. (*) P < 0.05 for wild-type lungs;( ns) not significant for 
Clock A19 lung; one-way ANOVA for the effect of time. [F] Temporal CLOCK binding to E-boxes in Nrf2 and Dbp gene promoters in 
wild-type and Clock A19 lungs using CLOCK-specific ChIR IgG served as a negative control. Data (mean ± SEM) were expressed as 
percent input. (*) P < 0.05; (ns) not significant, t-test. (White bar) Light phase,- (black bar) dark phase. 



contrast, the mRNA rhythm of Nrf2 was abolished and 
remained low across all time points in the lungs of 
Clock A19 mice, consistent with a loss of rhythms in the 
clock gene expression (Fig. 2E ; Supplemental Fig. SI). 
Consistently, knockdown of Bmall by siRNA signifi- 
cantly reduced the levels of Nrf2 mRNA in Rati fibro- 
blasts (P < 0.05) (Supplemental Fig. S2). In contrast, in 
synchronized MEFs from arrhythmic Ciyl~^~/Cry2~^~ 
mice, a model in which the negative arm of the clock is 
disrupted (Supplemental Fig. S3A,B ; van der Horst et al. 
1999), we observed constitutively high protein expression 



levels of NRF2 as well as mRNA levels of its two target 
genes (Gclm and Gsta3) (Supplemental Fig. S3C,D). 

To verify the function of the putative E-box element in 
the Nrf2 promoter in vivo, temporal chromatin immuno- 
precipitation (ChIP) assays were performed in the wild- 
type lungs using a CLOCK-specific antibody (Yoshitane 
et al. 2009). As a positive control, we first observed a 
time-dependent CLOCK binding to the E-box region of 
the Dbp promoter, a known CLOCK target (ZT10 > ZT22 ; 
P < 0.05) (Fig. 2F), as reported previously in the liver 
(Ripperger and Schibler 2006). Rhythmic CLOCK occu- 
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pancy was also observed in the E-box region of the Nrf2 
promoter (Fig. 2F) but not in the distal region of the 
promoter (Supplemental Fig. S4A). There was a signifi- 
cantly higher binding of CLOCK at ZT4 as compared 
with ZT16, coinciding with the peak and trough expres- 
sion of Nrf2 mRNA, respectively. Moreover, CLOCK 
binding to the Nrf2 promoter was arrhythmic and re- 
mained at an intermediate level in the Clock A19 lungs 
(Fig. 2F). Together, these results demonstrate that the 
Nrf2 gene is directly regulated by the core clock compo- 
nents in vitro and in vivo through the conserved E-box 
element in its promoter. 

Rhythmic promoter binding of NRF2 drives 
oscillations of ARE-regulated antioxidant genes 
involved in glutathione homeostasis 

Upon nuclear translocation, NRF2 protein binds to the 
AREs in the promoters of its target genes and thus 
maintains both basal and inducible expression of many 
antioxidant genes (Nguyen et al. 2009; Malhotra et al. 
2010). The rhythmic nuclear NRF2 levels prompted us to 
examine whether the expression of ARE-regulated NRF2 



target genes is also rhythmic. Temporal quantitative PCR 
(qPCR) of the mouse lungs demonstrated robust endoge- 
nous circadian rhythms for a number of NRF2 target 
genes, including Gclm, Gsta3, and Hmoxl (P < 0.05) (Fig. 
3 A; Supplemental Fig. S5), as well as a core clock gene, 
Bmall (P < 0.001) (Supplemental Fig.S5). To investigate 
the hypothesis that NRF2 may be rhythmically recruited 
to the AREs in the promoters of antioxidant genes in vivo, 
we performed temporal ChIP assays in the mouse lungs 
using an anti-NRF2 antibody. Here, we detected strongly 
rhythmic NRF2 binding specifically to the promoter 
regions spanning the well-characterized AREs (Malhotra 
et al. 2010) in Gclm and Gsta3 gene promoters (Fig. 3B ; 
Supplemental Fig. S4B). We observed significantly higher 
promoter occupancy at ZT0 and ZT6 (light phase) and 
lower occupancy at ZT12 and ZT18 (dark phase). The 
peak of NRF2 binding during light phase is in agreement 
with its maximal nuclear expression and the peaks of 
Gclm and Hmoxl mRNAs in the lung. Interestingly, 
despite maximal NRF2 binding to the Gsta3 promoter 
during the light phase, Gsta3 mRNA peaked during early 
subjective night, suggesting that additional phase delay 
mechanisms may be at play, such as the availability of 
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Figure 3. Circadian rhythm of NRF2 binding to ARE drives rhythmic oscillations of genes involved in glutathione synthesis 
and utilization. (A) Temporal mRNA levels of NRF2 targets in wild- type (WT) mouse lungs (DD). Data (mean ± SEM) were normalized 
to Gapdh, and the lowest expression was set as 1. [*)P< 0.05, one-way ANOVA. [B) Rhythmic NRF2 occupancy on AREs of antioxidant 
gene promoters in wild-type mouse lungs by NRF2-specific ChIP. The position of each ARE in relation to the TSS is shown. Mean ± 
SEM. (*) P < 0.05; (**) P < 0.01, t-test. (White bars) Light phase; (black bars) dark phase. (C) Temporal mRNA profiles of NRF2 target 
genes and clock gene Dbp in MEFs from wild-type and Nrf2 knockout (KO) mice following serum shock. Harvest of cells started at 28 h 
after serum shock. Data (mean ± SEM) were normalized to Gapdh. (*) P < 0.05 for all six genes in wild type; not significant for Nrf2 
knockout MEFs (except Hmoxl and Dbp; [*] P < 0.05), one-way ANOVA for the effect of time. (D) Temporal levels of reduced 
glutathione (GSH) in MEFs from wild-type and Nrf2 knockout mice following serum shock. Data (mean ± SEM) were normalized to 
cellular counts and quantified using standard curve for reduced GSH. (**) P < 0.01; (ns) not significant, t-test. 
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cof actors, as shown for the transcription factor BMAL1 
(Rey et al. 2011), or post-transcriptional processing 
(Patel et al. 2014). The difference in the timing of 
maximal expression may be necessary for the specific 
function of individual NRF2 targets. 

To investigate whether NRF2 is required for the rhyth- 
mic expression of antioxidant genes, we performed tem- 
poral qPCR for a panel of known NRF2 target genes 
in clock-synchronized MEFs from wild-type and Nrf2 
knockout mice. Here, we observed significantly reduced 
basal expression of many antioxidant genes in Nrf2 
knockout MEFs, as reported previously (Malhotra et al. 
2010). Interestingly, the rhythmic patterns of genes 
coding for enzymes involved in glutathione synthesis 
(Gclm and Gclc) and glutathione utilization (Gsta3 
and Gsta4) were either severely dampened or abolished 
(Fig. 3C). In contrast, Hmoxl, an inducible (but not basal) 
NRF2 target, remained robustly rhythmic in Nrf2 knock- 
out MEFs. Consistent with these results, the level of total 
glutathione (GSH) showed a significant temporal varia- 
tion in wild-type MEFs but remained constitutively low 
in Nif2 knockout MEFs (Fig. 3D). Importantly, the rhyth- 
mic expression of Dbp in Nrf2 knockout MEFs (Fig. 3C) 
and the behavioral free-running rhythms appear pre- 
served in Nrf2 knockout mice (Supplemental Fig. S6). 
These findings support a direct NRF2-dependent rhyth- 
mic control of its downstream targets rather than due to 
feedback actions of NRF2 on the molecular circadian 
clock. All together, these data demonstrate that the periodic 
NRF2 activity regulates the rhythmic GSH homeostasis 
through temporally controlling a key set of genes involved 
in glutathione synthesis and utilization. 

Time-of- day- dependent lung response to bleomycin 
challenge is coupled to temporal NRF2 activity 

The NRF2 antioxidant pathway is known to play a critical 
role in combating fibrotic injury (Cho et al. 2004; Walters 
and Kleeberger 2008). We therefore hypothesized that the 
rhythmic NRF2 activity may gate the tissue-protective 
responses to a fibrotic challenge in a time-of-day-depen- 
dent manner. To test this, we used an established in vivo 
mouse model of bleomycin-induced lung fibrosis (Walters 
and Kleeberger 2008). Bleomycin was administered intra- 
tracheally close to either the peak or nadir (ZT0 or ZT12, 
respectively) of NRF2 protein expression (see Fig. 1 A-C). 
As expected, 7 d after bleomycin challenge, characteristic 
histological changes, including areas of inflammatory cell 
infiltration and alveolar wall thickening (data not shown), 
were observed at both time points. In addition, we also 
observed a significant increase in collagen accumulation 
around the peri-bronchial regions, as revealed by Masson's 
Trichrome staining of lung sections and quantification 
of fibrosis using the Ashcroft scoring system (Fig. 4A ; 
Ashcroft et al. 1988). However, the degree of lung fibrosis 
was clearly dependent on the time of treatment and in- 
versely correlated with NRF2 protein levels. A signifi- 
cantly higher fibrotic score was observed when bleomycin 
was delivered at ZT12 as compared with ZT0 (Fig. 4B), 
associated with higher induction of fibrotic gene markers 



(Timpl, Colla2, and Mmp3) (Supplemental Fig. S7). Bleo- 
mycin challenge is known to activate NRF2-dependent 
antioxidant tissue responses (Cho et al. 2004). Following 
an acute bleomycin challenge (4 h), the induction of 
mRNAs for Gclc, Gsta3, and Hmoxl showed a clear 
time-of-day dependence, with a significantly higher in- 
duction at ZT0 than at ZT12 [P < 0.01 for all three genes) 
(Fig. 4C), consistent with the hypothesis of a clock 
"gated" antioxidant tissue response. 

To further establish the role of NRF2 rhythms in time- 
of-day susceptibility to bleomycin-induced fibrosis, we 
investigated whether the timed administration of sulfo- 
raphane (SFN ; a direct NRF2 activator) can reduce the 
degree of bleomycin-induced fibrosis seen at the circadian 
nadir in NRF2 levels. To this end, wild-type mice were 
pretreated with either vehicle or SFN administered in- 
traperitoneally at ZT6 (and thereafter at the same time 
point on days 1, 3, and 5) prior to a single challenge with 
bleomycin at ZT12 and assessed 7 d later. As a control, 
mice treated with SFN alone (at ZT6) showed increased 
NRF2 protein levels in the lungs at ZT12 as compared 
with vehicle treatment (data not shown). Here, the degree 
of bleomycin-induced lung fibrosis was significantly re- 
duced by the timed SFN treatments [P < 0.05) (Fig. 4D). 
Consistent with these findings, the bleomycin-induced 
expression of fibrotic markers (Colla2, Timpl, Mmp3, 
and Ctgf) was significantly blunted by SFN treatment 
(Fig. 4E). All together, these results strongly support a 
role for NRF2 rhythms in coupling antioxidant response 
to time-of-day susceptibility to bleomycin-induced lung 
fibrosis. Moreover, these findings suggest that timed 
SFN treatment might protect the lungs from oxidative 
challenges, especially those that fall within the nadir of 
NRF2 activity. 

Altered NRF2 /GSH pathway activity in the Clock A19 
lungs is associated with increased markers 
of oxidative damage and fibrotic injury 

Our data indicate that rhythmic NRF2 activity is coupled 
to the clock gated protective lung responses. To further 
establish the requirement of the functional circadian 
clock, we used the Clock 119 mouse model to investigate 
changes in the NRF2 pathway activity. Consistent with 
the loss of Nrf2 mRNA rhythmicity in the Clock A19 
lungs, NRF2 protein expression was constitutively low 
and arrhythmic (Fig. 5A). Moreover, ChIP assays revealed 
an arrhythmic and low NRF2 occupancy in the ARE 
region of the promoters for Gclm and Gsta3 (Fig. 5B). 
While qPCR revealed robust time-of-day-dependent mRNA 
expression of Gclm and Gsta3 in wild-type lungs, their 
expression in the Clock A19 lungs remained constitutively 
lower, at levels similar to or below the trough expression 
in wild- type lungs (Fig. 5C). 

In an attempt to test the hypothesis that the Clock A19 
lungs were more prone to bleomycin-induced fibrotic 
injury, unexpectedly, we identified a spontaneous fibrotic - 
like phenotype in Clock A1 9 lungs even without a bleomycin 
challenge. We found a significantly increased deposition of 
collagen fibers around the bronchioles in Clock A19 com- 
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Figure 4. Time-of-day-dependent lung response to bleomycin challenge is coupled to temporal NRF2 activity. (A) Representative 
mouse lung sections were stained with Masson's Trichrome solution to reveal collagen fiber deposition (blue staining) and overall 
fibrosis. Wild-type (WT) mice were treated with bleomycin [n = 8) or saline [n = 5) at either ZTO or ZT12 and, 7 d later, were sacrificed 
for tissue processing. Bar, 200 fxM. [B) Lung sections in A were scored for fibrosis according to the Ashcroft scoring system. Mean ± 
SEM; (***) P < 0.001, saline versus bleomycin, t-test; (*) P < 0.05, ZTO versus ZT12, two-way ANOVA. (C) Fold induction in wild-type 
lung mRNA levels of NRF2 targets [Gclc, Gsta3, and Hmoxl) following bleomycin delivered at ZTO or ZT12. Data (mean ± SEM) were 
expressed relative to saline controls at the respective time points. (**) P < 0.01, two-way ANOVA. [D,E) SFN treatment partially rescues the 
lung fibrosis phenotype, as revealed by Masson's Trichrome staining (with lower Ashcroft score) and reduced gene induction for fibrotic 
markers [Timpl, Colla2, Mmp3, and Ctgf). Wild-type mice were treated with either saline [n = 3), bleomycin [n = 4), or SFN plus bleomycin 
(n = 5) and, 7 d later, were sacrificed for histological and gene expression analysis. Mean ± SEM; (*) P < 0.05; (**) P < 0.01, £-test. 



pared with wild-type lungs (P < 0.001) (Fig. 6A). Consis- 
tently, the expression of fibrotic marker genes [Mmp3 and 
Ctgf) was significantly elevated in Clock A19 lungs (Fig. 
6B). A reduced GSH level is known to be involved in the 
fibrotic injury and remodeling processes (Liu and Gaston 
Pravia 2009). To explore the underlying mechanisms for 
the observed lung phenotype, we compared the temporal 
levels of reduced glutathione, a major antioxidant output 
of the NRF2 pathway. In wild-type lungs, we found a 
significantly higher level of reduced GSH at ZTO com- 
pared with ZT12 (Fig. 6C). In contrast, Clock A19 lungs 
showed constitutively lower levels of reduced GSH, 
equivalent to the nadir level in wild- type lungs. Similar 
alterations in GSH level were also found in synchronized 
Clock A19 MEFs, which could partially be restored by SFN 
treatment (Supplemental Fig. S8; Higgins et al. 2009), sup- 
porting the involvement of NRF2. To determine whether 
the low GSH levels in Clock A19 lungs correlated with an 



increased oxidative burden, we temporally examined the 
lung protein carbonylation levels (a marker of oxidative 
damage). We observed a time-dependent change in the 
protein carbonyl levels in wild-type lungs, which were 
anti-phasic to the reduced GSH levels (Fig. 6D). In con- 
trast, protein carbonylation levels remained high at both 
time points in Clock A19 lungs, consistent with a consti- 
tutive decline in GSH levels. These data support the 
hypothesis that the functional clock is required in co- 
ordinating the temporal activation of the NRF2/GSH 
pathway in the lung, disruption of which contributes to 
the increased susceptibility of the lung to cyclic oxidative 
stress and profibrotic challenges. 

Discussion 

It has recently been suggested that circadian clocks may 
mediate the protective responses to oxidative stress 
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Figure 5. Loss of rhythmic NRF2 pathway activity in Clock A19 lungs. (A) Representative Western blots of diurnal NRF2 levels in total 
lung lysates from wild-type (WT) and Clock A19 mice. NRF2 densitometry (mean ± SEM) was normalized to tubulin and expressed 
relative to wild-type ZTO. (*) P < 0.05 for wild type; not significant for Clock A19 , one-way ANOVA for the effect of time. (White bar) 
Light phase; (black bar) dark phase. [B] Temporal (ZTO vs. ZT12) ChIP assays for the ARE regions of NRF2 targets [Gclm and Gsta3) in 
wild- type and Clock A19 lungs using NRF2-specific antibody. Data (mean ± SEM) were expressed as percent input. (*) P < 0.05; (ns) not 
significant, t-test. (C) Temporal mRNA levels of Gclm and Gsta3 in wild-type and Clock A19 lungs at ZTO versus ZT12. Data (mean ± 
SEM) were normalized to Gapdh. (*) P < 0.05; (ns) not significant, t-test. 



(Stangherlin and Reddy 2013; Patel et al. 2014). However, 
the underlying clock-controlled molecular mechanisms 
and their role in tissue-specific physiology are still poorly 
understood. In this study, we identified the circadian 
clock as an endogenous molecular regulator of the 
NRF2/glutathione-mediated antioxidant defense path- 
way in the mouse lung, which couples the protective 
antioxidant response to the time-of-day variation in 
susceptibility to oxidant injury and pulmonary fibrosis. 
Based on our findings, we propose a working model de- 
picting how the circadian clock exerts its regulation of 
protective antioxidant responses in the lung through the 
control of transcription factor NRF2 under physiological 
conditions (Fig. 7). The circadian transcription factors 
CLOCK and BMAL1 positively regulate Nrf2 transcrip- 
tion through temporal control of the E-box element in the 
Nrf2 gene promoter. NRF2 protein thus accumulates in 
a circadian manner and drives circadian transcription 
of a key set of antioxidant genes involved in glutathione 
metabolism (such as Gclm and Gsta3). The periodic 
NRF2 activity is exerted through its rhythmic recruit- 
ment to the AREs in the promoters of target antioxidant 
genes and mediates a time-coordinated protection against 
oxidative tissue injury and fibrotic damage. 

Indeed, using an in vivo bleomycin-induced pulmonary 
fibrosis model, we show for the first time that, in the 
lung, the acute induction of the NRF2-mediated antiox- 
idant response is temporally "gated" by the clock and 
coupled to a time-dependent resistance to fibrotic injury. 
Importantly, a scheduled pharmacological activation of 



NRF2 in vivo was able to significantly mitigate the 
decrease in resistance to fibrotic effects of bleomycin 
delivered at the circadian nadir in NRF2 levels. More- 
over, in the lungs of Clock A19 circadian mutant mice, the 
altered temporal regulation of NRF2 antioxidant activity 
was associated with diminished glutathione levels, in- 
creased oxidative protein damage, and a spontaneous 
fibrotic-like phenotype. All together, our findings high- 
light an important role of the circadian clock in the 
regulation of the NRF2/GSH-mediated antioxidant re- 
sponse and its coupling to oxidative/fibrotic injury in 
the lung. 

Our study builds on the recent reports implicating the 
role of the circadian clock in the regulation of tissue 
physiology through redox-dependent mechanisms 
(Kondratov et al. 2006; Lai et al. 2012; Wang et al. 
2012; Musiek et al. 2013; Peek et al. 2013). For instance, 
direct circadian regulation of antioxidant genes [Nqol 
and Aldh2) has been reported in the cerebral cortex, and 
the circadian regulation of NAD + -dependent histone 
deacetylase SIRT3 has been implicated in mitochondrial 
oxidative metabolism. Our findings reveal a novel aspect 
of circadian clock function in modulating pulmonary 
fibrotic responses, which is mediated through its control 
of the NRF2 transcription factor. 

Current research concerning NRF2 regulation has 
largely focused on the post-translational mechanisms, 
especially via a KEAP1 -based protein turnover in re- 
sponse to various electrophilic/oxidative stresses. In this 
model, NRF2 proteins are tethered to the cytoplasm by 
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Figure 6. A fibrotic-like injury in the lungs of Clock A19 mice is associated with increased oxidative damage. (A) Representative wild- 
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KEAP1 and targeted for proteasomal degradation (Hayes 
and McMahon 2009). It was therefore thought that only a 
small pool of NRF2 protein can escape from its degrada- 
tion in the cytosol and translocate to the nucleus under 
basal conditions. In response to oxidative stress, however, 
oxidative inactivation of KEAP1 and de novo synthesis of 
NRF2 drive sustained nuclear NRF2 accumulation (Itoh 
et al. 2010). Our data reveal a rhythmic nuclear abun- 
dance of endogenous NRF2 protein in tissues and cells 
that accompanies the rhythmic levels in the total pool of 
NRF2 protein. These results are supported by previous 
studies in the mouse liver and pancreatic (3 cells that 
identified Nrf2 as one of the clock-controlled targets (Xu 
et al. 2012; Lee et al. 2013). It is therefore conceivable that 
during the circadian peak of NRF2 protein, increased 
levels of NRF2 may saturate the binding capacity of KEAP1, 
thus allowing NRF2 to escape from the KEAP1 -mediated 
degradation and translocate to the nucleus, where it in- 
duces higher expression of antioxidant genes. In contrast, 



during the nadir phase of NRF2 expression, KEAP1 may 
exert tighter control of NRF2 activity through increased 
NRF2 degradation. In this way, the regulation of NRF2 
nuclear activity appears to be restricted to a particular cir- 
cadian window, gated by both the clock-controlled mRNA 
synthesis and the KEAP1 -based protein degradation. 

Previous ChlP-seq (ChIP combined with deep sequenc- 
ing) and microarray studies have revealed that, in addi- 
tion to stress-induced antioxidant gene expression, NRF2 
directly controls basal expression of numerous antioxi- 
dant genes via AREs in their promoters (Reddy et al. 2007; 
Malhotra et al. 2010). Here we show that, under physio- 
logical conditions in the lung, the protein levels and 
recruitment of NRF2 to gene-specific antioxidant gene 
promoters is under rhythmic clock regulation. Nr/2-de- 
ficient MEFs are not capable of maintaining circadian 
rhythms of a number of key antioxidant genes involved in 
glutathione homeostasis and show arrhythmic and re- 
duced levels of its major antioxidant output glutathi- 
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one. Moreover, in Clock A19 circadian mutant mice, the 
altered regulation of NRF2 activity was also associated 
with reduced expression of genes involved in glutathi- 
one metabolism as well as total glutathione levels. 
Since both the behavioral free-running rhythms in Nrf2 
knockout mice and Dbp mRNA rhythm in Nrf2 knock- 
out MEFs appear preserved, this highlights the essen- 
tial role of NRF2 in mediating rhythmic glutathione 
homeostasis. 

Changes in the lung NRF2 expression and GSH metab- 
olism as well as elevated protein carbonyl content are 
recognized as central features of many human chronic 
lung diseases, including idiopathic pulmonary fibrosis, 
cystic fibrosis, acute respiratory distress syndrome, and 
COPD (Rahman and MacNee 2000; Dalle-Donne et al. 
2003; Rangasamy et al. 2004; Reddy et al. 2007; Liu and 
Gaston Pravia 2009). Nrf2 knockout mice show an in- 
creased propensity toward stress-induced tissue patholo- 
gies, including bleomycin-induced pulmonary fibrosis 
(Chan and Kan 1999; Cho et al. 2004). We observed 
temporal variations in protein carbonylation in wild-type 
mouse lungs, which were inversely correlated with 
changes in the reduced GSH levels. We show a time-of- 
day-dependent variation in the susceptibility to bleomycin- 
induced lung fibrosis in vivo, which was inversely corre- 
lated with the daily rhythm in the NRF2/GSH-dependent 
antioxidant response. Moreover, this time-dependent di- 
urnal variation in sensitivity could be mitigated by the 
scheduled administration of the NRF2 activator SFN, 
highlighting the NRF2 pathway as a potential chronother- 
apeutic target to prevent or slow down the progression of 
oxidant-induced injury and fibrosis. In the lungs of 
Clock A19 circadian mutant mice, the altered temporal 
NRF2 activity was associated with constitutively lower 
levels of GSH, elevated protein oxidative damage, and a 
spontaneous fibrotic-like phenotype. We therefore specu- 
late that decreased glutathione levels and GSH-dependent 
detoxification contribute to the increased oxidative/ 
fibrotic damage in Clock A19 circadian mutant mice. Con- 



sistent with our findings, Bmall knockout mice showed 
accelerated aging associated with constitutively high tis- 
sue ROS levels and oxidative stress-induced senescence 
(Kondratov et al. 2006, 2009), which could be rescued using 
glutathione precursor NAC (Khapre et al. 2011). 

In conclusion, we identified the molecular circadian 
clock as an endogenous regulatory mechanism control- 
ling the rhythmic activity of the redox-sensitive tran- 
scription factor NRF2 in the mouse lung. We also reveal 
the physiological significance of this rhythmic control in 
modulating tissue susceptibility to oxidative injury and 
pulmonary fibrosis. Given the role of NRF2 in human 
fibrotic diseases, we envisage a scenario in which circa- 
dian misalignments (e.g., caused by genetic variation or 
aging) may compromise the rhythmic activation of the 
NRF2-mediated antioxidant defense, leading to increased 
oxidative damage. Moreover, once fibrosis develops, there 
is no way to revert the phenotype. Current therapy that 
aims to slow down the progression of fibrotic diseases is 
often ineffective (Liu and Gaston Pravia 2009); our study 
implicates the need for considering circadian timing 
mechanisms, including timed drug administration (chro- 
nopharmacology) (Levi and Schibler 2007). Therefore, 
a major challenge in the future will be to understand 
how these circadian mechanisms affect tissue function in 
human disease. 

Materials and methods 

Animals and tissue collection 

All experiments were conducted under the aegis of the 1986 
Home Office Animal Procedures Act (UK). Mice were main- 
tained on a standard maintenance chow under a 12-h light/12-h 
dark (12:12 LD) regimen. PER2::Luc and Clock A19 mice on a 
C57BL/6J background were generated by Professor Joseph 
Takahashi. Clock A19 mice were subsequently bred with PER2— Luc 
mice. For circadian tissue collections or experiments with 
bleomycin-induced fibrosis, 8-wk-old female C57BL/6J mice 
(Harlan Laboratories) were placed under 12:12 LD cycles for 
2 wk before their release into DD. Animals were sacrificed by 
cervical dislocation in complete darkness using an infrared 
viewer, and lung tissues were harvested at 4-h intervals, begin- 
ning at 39 h after the start of DD. All tissues (circadian and 
diurnal collections) were either freshly used or snap-frozen in 
liquid nitrogen and kept at -80°C until use. 

Reagents and plasmids 

Bleomycin was purchased from ENZO. SFN was purchased from 
LKT Laboratories. All other chemicals were purchased from 
Sigma unless specified otherwise. mCLOCK and hBMALl ex- 
pression plasmids and pGL3-mPerl : : luc were kind gifts from Dr. 
Kazuhiro Yagita. The lentiviral mPer2::\uc construct was pre- 
viously described (Lu et al. 2010). The generation of wild- type 
and AE-box miYr/2::luc constructs is described in the Supple- 
mental Material. 



Real-time bioluminescence recordings 

Bioluminescence from the organotypic lung slices of wild-type 
and ClockAl9 mice (on a PER2::Luc background) (Yoo et al. 
2004) or from MEFs with lentivirally transduced Per2: :luc 
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reporters was recorded by photomultiplier tubes (PMTs,- 
Hamamatsu) housed in a light-tight incubator at 37°C and 0% 
C0 2 as described previously (Meng et al. 2008). The lung tissue 
slices were cultured in the recording medium (phenol red-free 
DMEM supplemented with 10% fetal bovine serum, 3.5 mg/mL 
glucose, 25 U/mL penicillin, 25 |xg/mL streptomycin, 0.1 mM 
luciferin, 10 mM HEPES-NaOH at pH 7) on Millicell cell culture 
inserts (Millipore). Cellular clocks were synchronized by a well- 
established serum shock protocol using 50% horse serum (for 
2 h). Data were presented as photon counts per minute. 

ChIP assays 

Freshly harvested lung tissues were homogenized and cross- 
linked in 1% (v/v) formaldehyde. Lung nuclear lysates were 
sonicated (—500 bp) using Bioruptor (Diagenode), and 5% of each 
sheared chromatin sample was saved as an input control. 
Following preclearing, lung chromatin was incubated overnight 
at 4°C with the following antibodies: 1 |xg/mL rabbit anti-NRF2 
(Santa Cruz Biotechnology), 4 |xg/mL mouse anti-CLOCK (CLSP3) 
(Yoshitane et al. 2009; in-house), and 1 |xg/mL control rabbit or 
mouse IgGs (Millipore). Immunoprecipitated DNA fragments 
were captured using protein G magnetic Dynabeads (Life Tech- 
nologies). After a series of washes, proteinase K digestion, and 
reverse cross-linking, DNA was eluted and cleaned up (Qiagen). 
Real-time qPCR was performed using TaqMan-based primer 
pairs/probes (Supplemental Table 1) to amplify the promoter 
regions spanning the conserved E-box or AREs. The amount of 
precipitated DNA in each sample was quantified according to 
the percent input method, in which Ct values of each sample were 
normalized to the adjusted Ct values of their respective inputs 
(corrected for dilution factor) and expressed as percent input. 

Tissue glutathione (GSH) assay 

Frozen lung tissues (100 mg) were homogenized in ice-cold 5% 
sulfosalicylic acid (Sigma) to minimize oxidation of reduced 
GSH. The diluted acid supernatants were used for GSH mea- 
surements, and the acid-precipitated protein was used for protein 
determination using a BCA protein assay kit (Pierce). Free 
reduced glutathione was determined using the luminescence- 
based GSH-Glow glutathione assay (Promega), which is based on 
the conversion of a luminogenic probe in the presence of reduced 
GSH in a reaction catalyzed by GST. Luminescence measure- 
ments were performed in triplicates for each biological replicate 
and processed in a single run to aid comparative quantification. 
Relative luminescence units (RLUs) were corrected for back- 
ground luminescence and converted to nanomoles per milligram 
of protein using a standard curve for reduced GSH. 

Protein carbonylation assay 

To measure protein carbonylation, lung protein lysates were 
processed using the Oxyblot protein oxidation detection kit 
(Millipore). Briefly, the carbonyl groups in the SDS-denatured 
proteins were derivatized to 2,3-dinitrophenyl hydrazone (DNP- 
hydrazone) by reaction with 2,4-dinitrophenylhydrazine (DNPH) 
followed by SDS-PAGE. Immunoblot detection of protein car- 
bonyls was performed using rabbit anti-DNP primary antibody 
and goat anti-rabbit HRP-conjugated secondary antibody. As 
a negative control, proteins were derivatized using a derivatiza- 
tion control solution. Protein carbonylation levels were ana- 
lyzed using National Institutes of Health ImageJ software by 
determining the values obtained for total protein density within 
each protein lane and were normalized to tubulin density for each 
sample. 



Bleomycin-induced lung fibrosis model 

All animal procedures were carried out in accordance with the 
Animals Scientific Procedures Act (1986). Female wild-type 
(C57/BL6), PER2::Luc, or clock A19 mice (aged 8-12 wk) were 
dosed intratracheally with 1 mg/kg clinical-grade bleomycin 
(30 [iL volume,- ENZO) or vehicle (saline) at either ZT0 (lights on) 
or ZT12 (lights off). Tissues were collected either 4 h later for 
RNA analysis of acute gene transcription or 7 d later (at ZT6) for 
histological assessment and long-term gene expression. The 
lungs were lavaged post-mortem via an incision in the trachea, 
with 1 mL of 10 mM EDTA and 1% BSA. For acute gene 
expression assays, fresh lung lobes were dissected and snap- 
frozen. For histological samples, 1 mL of 4% PFA was delivered 
via the trachea to inflate the lungs. The trachea was clamped off, 
and the lungs were removed en bloc and fixed in 4% PFA 
overnight at 4°C. After fixing, the middle right lobe was removed 
for RNA extraction, while the left lobe was processed for 
histological assessment. Lung slices at 5 |xm were stained with 
Masson's Trichrome solution following the standard procedure 
within our core histology facilities (University of Manchester). 
The histological slides were imaged and scored numerically 
using the Ashcroft scoring system by two experienced re- 
searchers blinded to the treatment groups. For rescue studies, 
10 mg/kg SFN (LKT Laboratories) or 10 mg/kg vehicle (saline) 
was administered intraperitoneally at ZT6 prior to a bleomycin 
challenge at ZT12 and then every other day for 7 d before tissue 
processing. 



Immunocytochemistry and immunohistochemistry 

Immunostaining was performed as previously described (Pekovic 
et al. 2011). Briefly, cells were fixed with 4% formaldehyde at 
12-h intervals after stimulation with 50% horse serum. For 
immunocytochemistry, rabbit anti-NRF2 antibody (sc-722, 
Santa Cruz Biotechnology) was applied at 1:200 overnight at 
4°C. The secondary antibody was donkey anti-rabbit IgG conju- 
gated to Alexa Fluor 555 (Molecular Probes). Nuclei were 
visualized by DAPI using VectaShield Hard-Set (Vector Labora- 
tories). Cells on cover slips were imaged using a Zeiss Axiovert 
40 CFL inverted microscope (Carl Zeiss). Intensity of NRF2 
immunostaining was quantified using ImageJ software by mea- 
suring the integrated density of fluorescence in 50-100 cells after 
correcting for background fluorescence, and the results were 
confirmed in three independent experiments. 

For immunohistochemistry, lung tissues were fixed with 4% 
formaldehyde, paraffin-embedded, and sectioned at 5 ixm. Slides 
were dewaxed, rehydrated, and processed for antigen retrieval. 
Endogenous peroxidases were blocked in 3% H 2 0 2 for 30 min. 
Blocking was carried out for 1 h in 2% normal goat serum, after 
which sections were incubated in anti-NRF2 antibody at 1:50 
dilution overnight at 4°C in a humidified chamber. Slides were 
incubated in biotinylated secondary donkey anti-rabbit antibody 
(Vector Laboratories) diluted 1:800 in 1 % BSA. ABC solution and 
the DAB (3'3'-diaminobenzidine) substrate kit (Vector Laborato- 
ries) were used to visualize positive NRF2 staining. To confirm 
the antibody specificity, immunohistochemistry was performed 
on control sections without the addition of either primary or 
secondary antibodies. 



Statistical analyses 

Data were evaluated using Student's t-test, one-way ANOVA 
with Tukey test, or two-way ANOVA for multiple comparisons 
as indicated. Results are presented as mean ± SEM from at least 
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three independent experiments. Differences were considered 
significant at the values of P < 0.05 (*), P < 0.01 (**), and P < 
0.001 (***). 
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